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It is known experimentally that at not very large filling factors v the quantum Hall conductivity 
peaks corresponding to the same Landau level number N and two different spin orientations are well 
separated. These peaks occur at half-integer filling factors v = 2N + 1/2 and v = 2N + 3/2 so that 
the distance between them hv is unity. As v increases hv shrinks. Near certain N = N c two peaks 
abruptly merge into a single peak at v = 2N + 1. We argue that this collapse of the spin-splitting at 
low magnetic fields is attributed to the disorder-induced destruction of the exchange enhancement 
of the electron (/-factor. We use the mean-field approach to show that in the limit of zero Zeeman 
energy 6v experiences a second-order phase transition as a function of the magnetic field. We give 
explicit expressions for N c in terms of a sample's parameters. For example, we predict that for 
high-mobility heterostructures N c = 0.9dn 5 ^ 6 n l 1 ^ 3 , where d is the spacer width, n is the density 
of the two-dimensional electron gas, and n t is the two-dimensional density of randomly situated 
remote donors. 

PACS numbers: 73.40.Hm 



I. INTRODUCTION 



A characteristic feature of the quantum Hall effect is 
the appearance of peaks in the diagonal conductivity 
<7 XX as magnetic field is varied. The conventional ar- 
gumentation for this is as follows. Consider the den- 
sity of states diagram for a two-dimensional electron gas 
(2DEG) in a perpendicular magnetic field (Fig. la). It 
consists of disorder broadened Landau level subbands 
(LLS): |,0 J,, 1 |,1 J., etc. Here arrows stand for the 
two spin orientations. When the magnetic field is de- 
creased, the Fermi level £f consequtively passes through 
these subbands. Whenever the Fermi energy coinsides 
with the center of some LLS, where there is a delocalized 
state, the peak in the dissipative conductivity occurs. 
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For the case shown at Fig. la where the width of LLS 
is smaller than the distance between them, peaks are po- 
sitioned at half-integer values of the average filling factor 
v: v = 2 A + 1 ± i, where A = 0, 1, . . . is the Landau 
level number, and the plus/minus sign corresponds to 
the spin down/spin up orientation. Experimentally, this 
periodic positioning of the peaks is observed at small v 
(high fields). In other words, at small v two peaks for 
the same A are well separated. The distance bv between 
them is approximately unity. At larger v (lower fields) 
the periodicity of the peak positions is violated. Peaks 
are found to appear in close pairs at v = 2A + 1 ± \bv 
with bv < 1. The experiment suggests that starting from 
a certain A = A c , this separation becomes so small that 
only single peaks instead of pairs are seen at odd integer 
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FIG. 1. The density of states diagram where the two upper LLS are shown 
as disorder-broadened peaks. The ratio of the shaded area to the total area 
of one LLS peak gives 6v. If the width of LLS is smaller then their separation 
As (a), then 6v is close to one. In the opposite case (b) it is much smaller. 
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Returning back to the density of states diagram, the 
smallness of bv can be interpreted as a result of the over- 
lapping of LLS (Fig. lb). The larger the overlap is the 
smaller the separation between the spin peaks become. 

It is worth mentioning that the possibility to observe 
a pair of peaks instead of a single one depends not only 
on their separation but also on the temperature because 
as the temperature increases, the conductivity between 
the peaks rapidly rises. At low temperatures it is, the 
most likely, due to the variable-range hopping. This 
mechanism was discussed by Polyakov and Shklovskii 3 . 
Recently, Polyakov and Raikh 4 studied the conductiv- 
ity between the peaks in a higher temperature range. In 
this paper we consider the case of zero temperature where 
peaks are very narrow and concentrate on the peak posi- 
tions only. However, we have to face the fact that all the 
experiments are performed at finite temperatures, and 
small bv can not be measured accurately. Nevertheless, 
we assume that the point bv = 1/2 can be reliably iden- 
tified and define N c as the smallest peak number where 
bv > 1/2. 

It appears from the experiment that bv is almost con- 
stant (= 1) away from N c but decreases rapidly in the 
very vicinity of N c . Usually, as N changes from, say, 
N c - 1 to N c + 1, bv drops from bv > 0.75 to bv < 0.25. 
For GaAs/Gai_ x Al x As heterostructures N c is usually in 
the range 1 — 12, and grows with the mobility and the 
setback distance of the doped layer 1,2 . 

In this paper we argue that the sharp change in bv, 
or in other words, the collapse of the spin-splitting at 
low magnetic fields, is attributed to the disorder-induced 
destruction of the exchange enhancement of the electron 
(/-factor 5,6 . In GaAs the contribution of the exchange in- 
teraction to the energy splitting between the LLS (the ex- 
change gap) exceeds the bare Zeeman energy Z = </oA«B-B 
by a factor of twenty 7,8 (here hb is the Bohr magneton, 
B is the magnetic field, and go is the bare electronic g- 
factor). This exchange gap can be defined as the energy 
necessary for creation of a well separated pair of an elec- 
tron in the spin down and a hole in the spin up LLS. 
The size of the exchange gap is the largest at odd in- 
teger filling factors where in the absence of disorder the 
spins at the upper LLS are completely polarized. We de- 
note this maximum value by Eo and the ratio Eo/Tiu> c 
by a. The parameter a is related to the (/-factor via 
a = i(g — go)(m/mo), where m and mo are the effec- 
tive and the bare electron masses, respectively. Exper- 
imentally 7,8 g ~ 6-7; therefore, a ~ 0.25. From the 
theoretical side, the calculation of the exchange gap for 
arbitrary Landau level number N was performed already 
in the early paper of Ando and Uemura 6 (see also the 
paper by Smith et al 9 and references therein). However, 
the explicit expression was given only recently by Aleiner 
and Glazman 10 who showed that a becomes independent 
of N when N gets large: 



a= ln(2^«B) ; JV>jbpaB>1) (l) 

irk^a^ 

where kp is the Fermi wave- vector and gib is the Bohr ra- 
dius in the semiconductor. In realistic situation k^a^ ~ 
1, so Eq. (1) may not apply literally. Nevertheless, for 
fcpciB = 1 Eq. (1) gives the value of a « 0.22. In our 
theoretical arguments we assume that k^a^ 1, but we 
will keep in mind that this parameter is close to unity 
when we make our estimates. 

So far we have not allowed for a possibility of the spin 
orientation being different from parallel or antiparallel 
to the magnetic field. It is known that by limiting our 
consideration in such a way we can reach incorrect con- 
clusions for v < 2. For example, Sondhi et al 11 showed 
that the lowest energy excitations for v = 1 are not 
the electron-hole but skyrmion-antiskyrmion pairs. (The 
skyrmion is a particular spin texture whose local spin 
orientation smoothly rotates from a point to point in 
space, the whole picture resembling a hedgehog. The 
antiskyrmion is a similar texture but with the opposite 
direction of rotation). Subsequently, however, Wu and 
Sondhi 12 found that for larger filling factors (v = 3, 5) the 
energy of these excitations 2Ask exceeds that of electron- 
hole pairs Eo- Therefore, already for these filling fac- 
tors spin textures are not the relevant physical states. 
Moreover, it can be shown that for large v the ratio 
2Ask/-E"o becomes much larger than unity 10 . For this 
reason throughout the paper where we are dealing with 
N > 1 we assume that the spin orientation can be only 
up or down. 

We argue that in the limit of zero Zeeman energy, in a 
mean-field approximation bv experiences a second-order 
phase transition as a function of the magnetic field. The 
possibility of such a transition has been mentioned ear- 
lier (Ando and Uemura 6 , Yarlagadda 13 , MacDonald and 
Yang 14 ). However, no explicit predictions for the transi- 
tion point N c were made. In order to make such predic- 
tions we have to choose realistic models for the disorder 
potential. Below we consider two limiting cases of short- 
range and long-range disorder. These cases supposedly 
describe low-mobility and high-mobility samples, respec- 
tively. 

In Sec. II we investigate the case of short-range disor- 
der and develop the mean-field description of the phase 
transition. The number N c in this case can be expressed 
in terms of the single-particle scattering time r: 

2a 2 n fi 

c = ~^h £FT ~ 10 10 cm- 2 10 6 cm 2 /Vs' ^ 

where in the last equation we used a = 0.25 and the 
fact that for the short-range disorder the single-particle 
scattering time r and the transport scattering time r tr , 
which enters the expression for the mobility /j, = er tr /m, 
are equal. 

Eq. (2) is in agreement with empirical observations 
that for samples with mobilities in the range jj, < 
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50,000cm 2 /Vs and the 2DEG densities n of order 2.0 • 
10 n cm _2 , it is usually only # = 0,1 peaks that are spin- 
split. 

In Sec. Ill and IV we consider a different mo- 
del, which applies to high-mobility modulation-doped 
GaAs/Gai_ x Al x As heterostructures. In this model the 
disorder potential is created by a plane of randomly po- 
sitioned ionized donors set back from the 2DEG by a dis- 
tance d. The Fourier harmonics of the disorder potential 
with wavelengths larger than d do not reach the 2DEG. 
Therefore, for large d the disorder potential is long-range. 
Its amplitude is determined by fluctuations bn^(r) of the 
two-dimensional donor density around its average value 
(«d)- Due to the Coulomb interactions between charged 
donors these fluctuations are usually much smaller than 
in the case of completely random distribution of donors. 
Following Refs. 15, 16 we assume that the fluctuations 
can be described by the correlation function 



( ( 5n D (T- 1 )(5n D (T-2)> = n^r-i - r 2 ), 



(3) 



where n 8 - <C («d) is the density of "uncorrelated" donors. 
This new parameter n 8 - depends on d and a "freeze-out" 
temperature 15-17 , and is not explicitly related to the ac- 
tual donor density. 

To find the disorder potential we also have to look at 
the screening properties of the 2DEG. To this end we 
has extended the theory of the non-linear screening (see 
Refs. 15, 18-21) to the case of weak magnetic fields (large 
N). A main new ingredient here was taking into account 
the screening properties of the lower N completely filled 
Landau levels. 

Using thus calculated disorder potential distribution, 
we then found the critical number N c . Depending on n, 
it is given by the following three expressions: 



STraV^n- 1 = 
n 1 / 2 In 2 (8wnal) 
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(4a) 
(4b) 
(4c) 



where is such a value of n that N c given Eq. (4a) 
matches the one given by Eq. (4b), i.e., is the solu- 
tion 22 of the transcendental equation 



n*/rii = [0.4/a(nH,)] a 



(5) 



Following Ref. 16 we have estimated n 8 - to be of the or 
der 0.1 1 • 10 n cm- 2 , so that 0.01 < n 8 a 2 < 0.1. Thii 
results into the values of n^a^ in the range 0.1 



0.3 or 



n* = 1.0-^3.0-10 cm . Having in mind that the exper- 
imental range of n is 0.4^-4 • 10 11 cm -2 , we conclude that 
Eq. (4b) is probably the most oftenly realized in the ex- 
periment. However, for large electron densities Eq. (4a) 
still may apply. 



Our predictions can be verified in detail on gated het- 
erostructures 23 . Changing the gate voltage, it is possible 
to vary n in a controlled manner and study the different 
regimes defined by Eqs. (4a-4c). Changing both the gate 
voltage and the magnetic field, one can fix N and study 
the collapse of a particular spin-split peak 24 . This may 
be a better way to study the critical behavior experimen- 
tally compared to the examining of bv at discreet values 
N . 

Concluding the introduction, we note that the collapse 
of the spin-splitting can be represented in terms of the 
global phase diagram for the quantum Hall effect origi- 
nally suggested by Kivelson et al 25 for the spinless case. 
The modified diagram is shown in Fig. 2. The purpose 
of this paper is to find the equation for the dashed line 
in this figure. 

The paper is organized as follows. In Sec. II we present 
the mean-field treatment of the problem and apply the 
results to the case of low-mobility samples. In Sec. Ill 
we study the case of high-mobility heterostructures with 
large electron densities: n/rii a~ 3 . We show that 
in this case the mean-field description is adequate. In 
Sec. IV we complement the mean-field picture by a mi- 
croscopic consideration and then derive from that the 
expression for N c in the case of moderate and small elec- 
tron densities n/rii <C a -3 . Finally, in Sec. V we make 
several concluding remarks. 
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FIG. 2. The global phase diagram for the quantum 
Hall effect as obtained in the mean-field approximation 
(schematically). The Zeeman energy is neglected. The 
dashed line is defined by the relation v = 2N C + 1 with N c 
given by Eqs. (4a-4c). 



II. MEAN-FIELD DESCRIPTION OF THE PHASE 
TRANSITION FOR THE SPIN-SPLITTING 

Our goal is to study the phase transition in the spin 
degree of freedom in the system of interacting electrons 
in weak magnetic field and in the presence of an external 
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disorder potential. As we will see this phase transition 
occurs when the amplitude of the disorder is still less 
than the cyclotron gap Tiu> c and therefore, only the two 
upper LLS are partially filled while all the lower 2N LLS 
are completely occupied. In this situation all the phe- 
nomena of interest take place at the upper Landau level. 
Therefore, a natural idea is to describe the electrons at 
the upper LLS by means of some effective Hamiltonian 
where the degrees of freedom of the lower LLS are in- 
tegrated out. This program was realized in Ref. 10. It 
uses the idea that the bare Coulomb interaction of the 
electrons at the upper LLS is strongly reduced due to 
the screening by lower LLS 26 . The fact that interaction 
is small suggests using perturbative methods, e.g, the 
Hartee-Fock approximation. 

Consider an electronic state at the upper LLS in the 
form of a Hartree-Fock function based on one-electron 
wave-functions proposed by Kivelson et at 27 (so-called 
coherent states). We will not need the explicit form of 
these wave-functions. It is sufficient to know that for 
large N the square of the absolute value (the probability 
density) of each wave-function of this type is not small 
only within a very narrow ring of width k^ 1 and radius 
R = v/k-p (the classical cyclotron radius). For a "ring" 
with the center at point ro the probability density at a 
point r is F(r — ro), where 



F(r) 



1 



2ttR 



6(r - R). 



(6) 



is the form-factor of the wave-functions. The ring cen- 
ters must be chosen such that their density be equal 
to «l = 1/(2 7I "'b)' where = \/v/kY is the magnetic 
length. Define the local filling factors v-\{r) and V[{r) 
as the fractions of occupied states with centers near the 
point r, and also vn = v^(r) + V[{r) as the total filling 
of the upper Landau level. 

Our effective Hamiltonian has the form 



H — Hi D t + -ffimp + Hz , 
-ffint = H ex + if H , 



(7) 
(8) 



where all the terms are expressed via the local filling fac- 
tors. For example, the exchange energy of the system 
ff™ is 



H ex = - X -nl jd 2 rd 2 r'G ex (r-r') 

x M^t<y) + n<>)n<y)]- 



(9) 



The kernel of this expression, G ex (r), is the energy of the 
exchange interaction of two "rings" whose centers are a 
distance r apart. The Hartree energy of the system Hn 
is given by 

Hh = \ni jd 2 rd 2 r'v N (r)G li (r - r')v N (r'), (10) 

where Gn(r) is now the energy of the direct interaction 
of two "rings" . The term ff; mp describes the interaction 
with random impurities: 



H 



imp 



n L Jd 2 rd 2 r'U(r')F(r' - r)v N {r), (11) 



where by U(r) we mean the impurity potential already 
screened by the lower completely filled LLS. The form- 
factor F here [see Eq. (6)] accounts for the fact that the 
energy of each one-particle state is determined by the po- 
tential averaged over the the cyclotron orbit 28 . To sim- 
plify the formulae we introduce an additional notation 
for this averaged potential: 



V(r) = d 2 r'U(r')F(r' -r), 



so that 



H 



imp 



n L d 2 rV(r)v N (r), 



Finally, there is the Zeeman term 

Hz = ^n L Z jd 2 r[v i {r) - v^r) 



(12) 



(13) 



(14) 



To complete the definition of the Hamiltonian we also 
have to provide the explicit expressions for G ex (r), Gn(r) 
and U (r). 

According to Ref. 10 



Gex(r) 



e 2 a B 



1 



2ttR r(l + r/£,)' 



< r < 2R, 



where 



-R. 



(15) 



(16) 



G ex (r) rapidly falls off beyond the distance £, which can 
be called the range of the exchange interaction. It is easy 
to verify that if N (A^ae) 2 , then £ ^> /b- Hence, for 
such N the number of electrons at the upper Landau 
level involved in the exchange interaction with the given 
one is large and we can indeed describe the distribution 
of electrons by continuous functions v-\{r) and V[{r). In 
the experimental situation k^a-B ~ 1, and therefore, we 
are referring to N ^> 1. 

The explicit formulae for Gn(r) and U(r) will not be 
used in the present Section. We postpone their discussion 
until Sec. III. 

The ground state distributions of v-\{r) and V[{r) can 
be inferred from the fact that H has the lowest value for 
a given total number of particles, which gives 



v^r) 
v i( r ) 



f[W(r) 
f[W(r) 



Z/2 
Z/2 



E^r)] 
Ei(r)] 



(17) 
(18) 



where 
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W(r) 



V(r) + n L Jd 2 r'G H (r 



r')v N (r') (19) 



is the potential, which includes the screening effects of 
both the lower completely filled LLS (the first term) and 
the upper partially filled ones (the second term), /(e) is 
the Fermi-Dirac distribution function, and is given by 



E\ 



dVG ex (r-r> T (r') 



(20) 



(and similarly for E^). Notice that for every position r 
of the ring center there are two distinct values of energy 
W(r)-Z/2-E ] {r) and W(r) + Z/2-E l (r) correspond- 
ing to the two orientations of the electronic spin. They 
define Landau level subbands. The energy separation 
As(r) between the subbands is a sum of two terms: 



A s = Z + E, 
E = E<i - E[. 



(21) 
(22) 



Eqs. (17,18) are in agreement with the density of states 
diagram (Fig. 1), which we discussed above. Using 
Eqs. (20,22) we can rewrite E in the form 



E(r) = d 2 r'G ex (r - r')[^(r') - v^r')}. 



(23) 



If both v-\{r) and V[{r) are uniform on the scale of £, 
then E is proportional to the difference in the subband 
populations 6 : 



E = Eoi^ - v{), 



(24) 



where Eq = Jd 2 rG ex (r) = ahu> c . Performing this inte- 
gration with the help of Eq. (15), we recover Eq. (I) 29 . 

Since we are interested in the conductivity peaks, from 
now on by v-\ and v^ in Eq. (24) we will understand the 
corresponding quantities for the position of the Fermi 
level at the center of N | LLS as shown in Fig. 1. Eq. (24) 
is one of the central equations in the mean-field theory 
of the spin-splitting phase transition. Another equation 
relates v-\ — i>[ to the separation bv between the conduc- 
tivity peaks: 



2 



(25) 



Indeed, it is easy to see that bv is the total area of the 
shaded strip (of width As) at Fig. 1, whereas v-\ — i>[ is 
only a half of this area. 

One more equation relates bv to As: 



bv 



p(e)de, 



(26) 



thus completing the system of the mean-field equations 
Eqs. (21,24,25,26). Here p{e) is a new quantity, which we 
call the form-factor of LLS. It has the dimensionality of 
inverse energy and has the direct relation to the density 



of states in the potential W(r). Namely, the contribu- 
tion to the to total density of states from, say, N | LLS 
centered at the energy — £j is given by ni,p{e + E^). 

Using Eqs. (21,24,25) we arrive at 



A= 



1 



E bv + Z. 



(27) 



The first (exchange) term is much larger than the sec- 
ond (Zeeman) one, provided bv is not too small. First, 
let us neglect the Zeeman energy. In this approximation 
the collapse of the spin-splitting can be seen the most 
vividly. We examine the effect of a finite Zeeman energy 
a bit later. 

We are going to solve the system of two equations (26) 
and (27). This procedure can be illustrated graphically 
(Fig. 3). We have a non-zero solution for bv only if 
Eq > /9 _1 (0). Therefore, the transition point is defined 

by 



E o = p- 1 (0). 



( a ) E <p(0) E >p(0) 



(28) 




(b) 



p(0)A s 

Magnetic field 
decreasing 




p(0)E 



FIG. 3. (a) hv is given by the intersection of the bold 
lines representing Eqs. (26) and (27). Non-zero solution 
exists only for Eo > /9 -1 (0). (b) The dependence of 6v 
on Eo has a typical second-order phase transition form. 
Non-zero Z smears the transition. To generate these 
plots we used the Gaussian form for the density of states: 
p(e) = p(0) exp( — 7T/9 2 (0)e 2 ). The solid and the dashed 
lines correspond to Z = and Z = 0.06i?o, respectively. 
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Near the transition bv can be estimated using Taylor 
expansion for p{e) in Eq. (26): 



bv : 



-24 


1/2 




1/2 


[p"(0)E 3 o \ 




[ E I 





(29) 



It is easy to see that our approach is very similar in spirit 
to such well-known mean-field theories as Stoner model of 
band ferromagnetism 30 or Weiss' theory of the molecular 
field. In the latter case the analogy is achieved via sub- 
stitutions "disorder" (p -1 ) — > temperature 31 , exchange 
energy (E) — ► magnetization. Following this analogy, 
the Zeeman term is mapped onto the external field. It is 
clear now that the effect of this term is the smearing of 
the transition (Fig. 3b). For instance, in the "paramag- 
netic phase" we have an analogue of Curie- Weiss law: 



bv 



2Z 



p-^O) -E ->Z. 



(30) 



Now let us address the question how reliable the mean- 
field approach is. It is based on the assumption that v-\ 
and V[ can be treated as spatially uniform. In reality, 
they fluctuate from point to point and in the immediate 
vicinity of the transition where these fluctuations exceed 
the average value of v-\ — v^, the mean-field theory def- 
initely breaks down. However, in this paper we do not 
wish to discuss the details of the critical behavior at small 
bv. Instead, we focus on not too small bv, say, bv > 1/2. 
The characteristic scale of the fluctuations in v-\ and v^ 
is set by the correlation length / cor of the disorder po- 
tential. As long as / cor <C ^, the fluctuations in v-\ — v^ 
are averaged out, and therefore, the mean-field theory is 
a very good approximation. For example, it is expected 
to work for low-mobility samples where the disorder is 
short-range, i.e., / cor <C /b (and we assume N (A^ae) 2 
so that /b "C £ as explained above). In this case we can 
express N c in terms of the sample mobility. Indeed, it 
has been shown 28 that for the short-range disorder the 
SCBA approximation is valid. In this approximation 



p(0) = h~ V(2/tt)(tK), 



(31) 



Together with Eq. (28) this leads to Eq. (2), which as 
we mentioned in the Sec. I, qualitatively agrees with the 
experimental data. 

In the next Section we discuss the case of high-mobility 
heterostructures with n/rii a~ 3 . As we will see, the 
simple theory presented in this Section gives the adequate 
description of the phase transition. 



III. HIGH-MOBILITY HETEROSTRUCTURES: 
LARGE ELECTRON DENSITIES 

In the preceding Section we applied our mean-field ap- 
proach to low-mobility samples where the disorder po- 
tential is short-range. Let us now turn to the case of 



high-mobility samples. In particular, we consider a het- 
erostructure described by a set of three parameters n, n;, 
and d as discussed in Sec. I. In this Section we will find 
the transition point N c assuming n/rii a -3 1. 

At first sight the problem is different from the one in 
the preceding Section because the disorder potential cre- 
ated by remote donors is now long-range with the corre- 
lation length of order d l B • This length also sets the 
scale for the variations of and v^ [Eqs. (17,18)]. How- 
ever, we will show below that the transition point is given 
by Eq. (4a); therefore, near the transition the range of 
the exchange interaction £ ~ 2iraR ~ a 3 dn/rii is much 
larger than d, and so on the scale of £ both v-\ and v^ are 
uniform and our mean-field treatment still applies. The 
problem now reduces to the calculation of the r.h.s of 
Eq. (28), i.e., to the calculation of the LLS width. This 
width is nothing else than the amplitude of W(r). More 
precisely, 



p(O)- 1 = y/2ir{W 2 ) 



(32) 



if W(r) has the normal distribution 28 . To calculate this 
quantity we need to know what the potential created by 
fluctuating donor density is, and how it is screened by 
the 2DEG. The former is characterized by exponentially 
suppressed short-scale harmonics (with the wave- vectors 
larger than d~ r ). The screening in its turn greatly re- 
duces the amplitudes of large-scale harmonics. As a re- 
sult, the harmonics with q ~ c? _1 dominate in the poten- 
tial, and so we can focus only on the screening at these 

Both the lower, completely filled LLS (the incompress- 
ible or i- liquid) and the upper, partially filled LLS (the 
compressible or c- liquid) participate in the screening, 
which is reflected in the fact that W(r) as given by 
Eq. (19) is a sum of two terms. The screening by the 
i-liquid is described by the dielectric constant 10,26 



£ = £s X 



l + R 2 q/a B , q<l/R 

l + 2/a B q, l/i?<g<A; F , 



(33a) 
(33b) 



where e s is the dielectric constant of the semiconductor. 
In the asymptotic expression for large q [Eq. (33b)] one 
recognizes the familiar zero magnetic field result 32 . This 
is natural because large q correspond to small distances 
</ _1 <C R at which the curvature of the cyclotron orbits 
may not be important. The dominant harmonics with 
q ~ g? -1 belong to this range because near the transi- 
tion £ d and, of course, R £, so that R ^> d. We 
will show that the contribution to the screening of the 
dominant harmonics from the c-liquid is negligible, and 
so 



W(q) « V(q), q ~ d~ 



(34) 



Therefore, for calculation of \J (W 2 ) we will need to con- 
sider only the screening by the lower LLS and only in the 
regime q ^> l/R described by Eq. (33b). Remind that 
V is the potential averaged over the cyclotron orbits and 
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in real space it is expressed in terms of the unaveraged 
potential U by means of Eq. (12). Now we rewrite this 
equation in (/-space: 

V(q) = F(q)U(q). (35) 

where F(q) can be deduced from Eq. (6) to be 

F(q) = Jo(qR), q < k F . (36) 

Due to the averaging over the cyclotron orbits, the am- 
plitude of V(r) is reduced compared to that of U(r) by 
a factor of order ^/d/R 28 . This is because every cy- 
clotron orbit samples approximately R/d uncorrelated 
areas of size d each. Another way to see that is to no- 
tice that f or th e dominant harmonics q ~ c? -1 we have 
\F(q)\ ~ ^/d/R. 

From Eqs. (32,35) we obtain 



dqq.J 2 (qR)S(q), 



(37) 



where S(q) is the structure factor of the disorder poten- 
tial defined via 



(U(qi)U(q 2 )) = (27r) 2 S(q)6( qi + q 2 ). 
Using this definition and Eq. (3) we find 

"r, 2 1/2 



so that 
1 



P(0) 



S(9) 



- qd 



2 1/2 

7r e z a B n,- 



R 



1/2 



Together with Eq. (28) this leads to 
N c = 2a 2 kpdn/rii, 



(38) 



(39) 



(40) 



(41) 



which is the same as Eq. (4a). 

Now we would like to justify our assumption that the 
screening by the c-liquid is negligible. We will show that 
this is true provided n/rii a~ 3 . 

The screening by the c-liquid may be both in the linear 
and the non-linear regimes depending on the amount of 
this liquid. This amount varies with the filling factor v, 
which gives rise to numerous screening regimes in differ- 
ent regions of the global phase diagram of Fig. 2. Since 
we are interested primarily in the positions of the con- 
ductivity peaks, we will focus here on a small part of this 
diagram near the phase boundaries defined by the equa- 
tion v = 2N C + 1 ± \bv{ni, n). As we move along such a 
boundary in the direction of the increasing disorder, we 
proceed from the point where the amplitude of W is zero 
and the spin-splitting is the largest, to the point where 
this amplitude is aTiu> c and the spin-splitting vanishes. 
Consequently, all the way along this path the amplitude 



of W never exceeds aTiu> c and moreover Tiu> c . This means 
that as long as N < N c , the c-liquid forms only at the 
two upper LLS (for the exception of rare places where 
W(r) is untypically large). 

To estimate the (averaged) potential created by the c- 
liquid we need to know the expression for Gn(r). It is 
given by 



Gh(t-) 



d 2 q 27re 2 
(2tt)2 qe(q) 



F 2 (q)e lq - r 



(42) 



where we took into account the form-factor of the wave- 
functions [Eq. (36)] and the screening by the 2N com- 
pletely filled lower LLS [Eq. (33)]. This integral was 
evaluated in Ref. 10 to be 



G H (r) 



e 2 a B 



2irRr 4ir 2 R 2 



3e 2 a B , ( R 
In — 



r / R 2 
B < r < R 



e 2 a B ^ (R 



(43) 



Among these three terms, the last one is just a constant 
and will have no effect in our consideration, and the sec- 
ond term is smaller than the first one. Retaining just the 
first term we discover a remarkable fact: the renormal- 
ized interaction is equivalent to the Coulomb law with an 
effective dielectric constant 



2ttR 
a B 



(44) 



Now the amplitude of the potential created by the c- 
liquid can be easily estimated. The amplitude of the 
variations in vn{v) does not exceed one; therefore, the 
amplitude of the potential is of the order of e 2 n\ J d/£* ~ 
e 2 a-Bn 3 / 2 d/N 2 , which is much smaller than \J (W 2 ) if 
n/rii ^> a -3 and N ~ N c [see Eqs. (40, 41]. Therefore, 
the screening is performed mainly by the i-liquid as we 
claimed above 33 . 

The expressions for N c obtained for the short-range 
disorder potential considered in Sec. II and in the present 
case appear to be different. Let us show that, in fact, they 
can be cast in a very similar form. Indeed, Eq. (37) can 
be rewritten in the form 



p(o)- 1 = h^fJJi 



(45) 



which differs from the short-range case [Eq. (32)] only by 
a numerical factor In this expression r is again 

the single-particle scattering time at B = 0: 



h 2 k 



F ^0 



dq,S(q). 



(46) 



With the help of Eqs. (28,45) N c can now be presented 
in the form almost identical to Eq. (2): 



N, 



-e F T. 



(47) 
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Therefore, the difference between the two considered 
cases is largely in the choice of measurable sample pa- 
rameters, through which we express N c . While for the 
short-range disorder potential we can extract r from the 
mobility, in the long-range disorder potential case we use 
n, rii, and d. 

Concluding this Section, we would like to mention that 
r can be estimated by analyzing the Shubnikov-de Haas 
(SdH) oscillations. Existing theories 32,34,35 predict that 
the deviation of the resistivity in a weak magnetic field 
from its zero field value should behave as 

Pxx " P ™ = -Ae-*'»° T x cos ( ^) . (48) 

There is a disagreement in literature about the value of 
the numerical coefficient A (see the discussion in Ref. 35) 
and problems with fitting the experimental data to the 
theoretical predictions 36 . Nevertheless, it is probably 
reasonable to assume that the peak number #sdH where 
the extrapolated to T = value of the r.h.s. of Eq. (48) 
becomes of order unity corresponds to lo c t ~ 7r, i.e., 
^VsdH ~ e^r/irk. Comparing with Eq. (47) we obtain 

N c ~ a 2 7r#sdH. (49) 

With a = 0.25 this gives N c ~ in reason- 

able agreement with the data by Coleridge et al. 2 where 
N c = 12 and #sdH ~ 40. The reason why there is such 
a simple relation between N c and #sdH is that the un- 
derlying physics is similar: the SdH oscillations develop 
when the width of the LLS becomes comparable with the 
energy separation Tiu> c of the Landau levels, and the col- 
lapse of spin-splitting occurs when this width is of the 
order of the spin subband separation akui c , which differs 
only by a numerical factor. 

In the following Section we are going to complement 
the mean-field description we employed so far by the con- 
sideration on the microscopic level. This will enable us 
to find N c for the case of moderate (and small) electronic 
densities n <C rii/a 3 . 



IV. HIGH-MOBILITY HETEROSTRUCTURES: 
MODERATE AND SMALL DENSITIES 

Consider now the case of moderate electron densities 
1 <C n/n l a~ 3 , where the predicted transition point is 
given by Eq. (4b). One can verify that near the transition 
we still have the inequality d <C R so that the averaging 
over the cyclotron orbits is still important for the calcu- 
lation of the macroscopic density of states 28 . However, 
the range of the exchange interaction £ ~ aR, is shorter 
than the lengthscale on which this density of states is 
formed (d). In this case the mean-field theory we used 
so far becomes insufficient, and there is a need for a new, 
local theory. 



A first step towards constructing such a theory will be 
to deepen our understanding of the already studied large 
density case. The starting point of our analysis was to 
associate the peaks in <r xx with the presence of the de- 
localized ("metallic") states near the Fermi energy. Now 
we have to somehow relate these metallic properties to 
local values of the filling factor v{r), which can be either 
integer (i-liquid) or non-integer (c-liquid). Clearly, the 
i-liquid where there is no gapless excitations, is unable 
to carry a dissipative current. It is "insulating" in this 
sense. Note that in the large density case the i-liquid oc- 
cupies almost the entire area of the sample. This is due to 
the fact that near the transition only the two upper LLS 
are typically partially filled and the potential created by 
the electrons at these LLS is small. The c-liquid exists 
only in the form of narrow channels between neighboring 
"islands" of the i-liquid. These channels are referred to 
as "bulk edge channels" . The peaks in <r xx appear when 
at certain values of the average v, the bulk edge channels 
form a percolating network, which can conduct the dis- 
sipative current. The transport properties of such a net- 
work have been considered by several authors 37,38 . For us 
it is important to identify the difference in the structure 
of this network for the "ferromagnetic" and "paramag- 
netic" phases. This difference is illustrated in Fig. 4a 
where we show the structure of the channel network for 
the periodic external potential 

^^^K^)™^)' (50) 

corresponding to the averaged potential V(r) of the same 
form: 

V(v) = U(r)J (R/d). (51) 

It follows from Eqs. (17,18) that in the "ferromagnetic" 
phase where As > 0, the total area occupied by the i- 
liquid is split into three regions where vn{v) takes values 
(both spin up and spin down LLS are empty), 1 (the 
spin up LLS is occupied and spin down LLS is empty), 
and 2 (both LLS are occupied). The bulk edge channel 
network consists of two subnetworks with different spin 
orientations: one is the channels with < vn < 1 (spin 
up) and the other is those with 1 < vn < 2 (spin down). 
The two subnetworks are disconnected as they are the 
two boundaries of region "1". The percolation through 
the two networks is achieved at different average v, and 
this is why there are two peaks in <r xx for a given N . 
Let us find when these peaks occur. Note that the spin 
up and spin down channels follow the contours of con- 
stant V(r) (or, more precisely, of constant W(r) but in 
this case they are the same). Only one such contour, 
namely, V = percolates through the sample (V(r) is 
symmetrically distributed around zero). Therefore, the 
peaks in <r xx correspond to £f = —E^ and £f = ~E[, 
which are the centers of the spin up and spin down LLS 
[see Eqs. (17,18)]. Thus, the network picture is consistent 
with the one we used previously, and gives an important 
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insight into the nature of the delocalized states at the 
centers of the LLS. 

On the contrary, in the "paramagnetic" phase where 
As = (Zeeman energy neglected) only the regions "0" 
and "2" are present (Fig. 4b). In this case the spin up 
and spin down networks coincide spatially, and the perco- 
lation through the two is achieved simultaneously, which 
means that bv is zero. 




FIG. 4. The distribution of v(r). Labels "0", "1", and 
"2" correspond to filling factors 2N , 2N+1, and 2N+2, re- 
spectively, (a) "Ferromagnetic phase" (6v > 0). The spin 
up bulk edge states (solid lines) are disconnected from the 
spin down ones (dashed lines). The space between them is 
region "1" (shaded), (b) "Paramagnetic phase" (6v = 0). 
Only the regions "0" and "2" are present and the two 
networks spatially coincide. 

In Sec. II and III we found that the transition is driven 
by the competition between the exchange interaction and 
the disorder. Now we can trace this competition on the 
microscopic level. The disorder is indiscriminative to- 
wards the electron spin. It tries to push the electrons, 
spin up and spin down all the same, away from the hills of 
the potential and pile them up in the valleys, thus tending 
to create doubly-empty ("0") and doubly-filled ("2") re- 
gions only. The exchange interaction, on the other hand, 
prefers the ferromagnetic spin ordering and leads to spin- 
flip processes at the spin down LLS. Since in the region 



"2" all the spin up states are occupied, these spin-flip 
processes occur via filling of the region "0" at the cost 
of an increase in energy due to the disorder potential. 
Thus, the region "1" is created, and the spin up and spin 
down bulk edge channels become spatially separated. 

In our example of the chess-board potential the perco- 
lating constant energy contour is the one, which goes 
through the saddle points. Therefore, the properties 
of the channel networks near the saddle points are ex- 
tremely important for the percolation (this is true for 
the random external potential as well). In fact, it is near 
a saddle point that the spin up and the spin down chan- 
nels split: one turns in the clockwise and the other in the 
counterclockwise direction (Fig. 4a). One can say that on 
the microscopic level the spin-splitting transition is the 
transition of the saddle points from a zero magnetization 
"0"-"2" state to a magnetized "0"-'T"-"2" state. Recall 
that bv differs from the average magnetization only by 
a factor of two [Eq. (25)]; therefore, bv is approximately 
given by the the ratio of the magnetized area at a sin- 
gle saddle point to the area of a chess-board cell. The 
former can be estimated as the square of the distance bx 
between the two regions "2" (or two regions "0") at the 
saddle point in the configuration shown in Fig. 4 (corre- 
sponding to v = 2N + 1). This leads to 

bv ~ (bx/d) 2 , (52) 

We see that for bv = 1/2 one should have bx of order d. 
We will use this fact later in this Section. 

Before we start to analyze the moderate electron den- 
sity case let us give another explanation why the mean- 
field description is expected to work in the large n case. 
The point is that in the "ferromagnetic" phase each mag- 
netized saddle point creates the "exchange field" , which 
supports the existence of its own magnetization and the 
magnetization of all neighboring saddle points. In the 
case of large n where the distance between the saddle 
points (d) is much smaller than the range of the inter- 
action (£), each magnetized saddle point is supported by 
the field of large number of other saddle points. There- 
fore, this is a collective phenomenon and the mean-field 
approach is a good one. 

Let us now resume the study of the moderate electron 
density case. In this case £ <C d and each saddle point has 
to support its magnetization individually. Another differ- 
ence from the large density case is that the potential cre- 
ated by the electrons at the upper LLS is not small, and 
the screening is now accomplished mainly by the upper 
LLS. In other words, both the Hartree and the exchange 
terms of the electron-electron interaction at the upper 
LLS are working against the disorder. The important 
role of the Hartree term was emphasized by Dempsey et 
a/ 39 who studied a similar "0"-"2" to "0"-'T"-"2" transi- 
tion for N = in an edge of a wide quantum wire. (See 
also the paper by Manolescu and Gerhardts 40 and refer- 
ences therein). The corresponding effect in our geometry 
is the "0"-"2" to "0"-"l"-"2" transition on the sides of 
the chess-board cells. We remind that our interest is 
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primarily in the saddle points. However, our treatment 
of the problem is similar in spirit; therefore, we briefly 
review the results obtained in Ref. 39. Let x be the co- 
ordinate in the direction transverse to the wire and the 
edges of the wire be located at x = and x = —b. We 
denote the potential providing the lateral confinement 
of the wire by V(x) (this choice reflects the similarity 
between this potential and the potential V(x,y) in our 
chess-board model). The evolution of the edge structure 
as V(x) is softened is as follows. For large confining fields 
dV/dx the edge is in "0"-"2" state (v = 2 in the interval 
— b < x < and v = everywhere else, see Fig. 5a). 
When the confining field gets smaller than the Hartree 
field given by 



dx 



(53) 



a) * t 
♦ ♦ 



b) 



x 



c) t 




x = 0, 



FIG. 5. The "0"-"2"to "0"-"l"-"2" transition in an 
edge of the wire 39 , (a) For a steep confining potential 
the edge is in "0"-"2" state, (b) As the confining poten- 
tial is softened, the region "1" appears and grows, (c) 
When the width of the region "1" reaches Ib , the c-liquid 
strips appear, (d) Eventually, the c-liquid strips become 
much wider than region "1". 



the region "1" arises spontaneously 39 and then grows 
(Fig. 5b). For a while the edge density profile remains 
step-like, i.e., the edge channels are infinitesimally nar- 
row. However, when the width of the region "1" becomes 
of order /b, the edge undergoes a second transition: the 
edge channels acquire a finite width (Fig. 5c). Using 
another language, the c-liquid appears. One can inter- 
pret this second transition as follows. It is the exchange 
interaction that stabilizes the step-like density profile 39 
because it provides the minimum of the exchange en- 
ergy H ex [Eq. (9)] for a given total magnetization. How- 
ever, this kind of profile does not simultaneously mini- 
mize the Hartree (i.e., electrostatic) part of the energy 
i?H [Eq. (10)]. Electrostatics tends to create a smoothly 
varying density, and starting from the point where the 
width of the region "1" gets comparable with the range 
of the exchange interaction, which is /b at the lowest 
Landau level, the exchange interaction is no longer able 
to counteract this tendency. With further softening of 
the confining potential (Fig. 5d), the c-liquid regions be- 
comes much wider than the region "I" 41 (this regime was 
considered by Chklovskii et a/ 42 ). 

In our chess-board model the situation is similar: when 
the interaction strength is comparable or larger than that 
of the disorder, the narrow bulk edge channels along the 
sides of the chess-board we were talking about earlier in 
this Section, become wide compressible regions occupy- 
ing a large portion of the total area. The saddle points 
undergo the same reconstruction. We are going to show 
below that for N close to N c where bv ~ 1/2, the size of 
the region "1" at the saddle points is much smaller than 
the size of the surrounding compressible areas (Fig. 6). 
In this case the profile of v^(t) in the cross-section going 
through the saddle point via labels "2"-"l"-"0" in Fig. 6a 
is similar to that shown in Fig. 5d. 



The fact that c-liquid regions are wide has qualita- 
tively new implications for the transport properties. The 
question now is whether or not this entire compressible 
area retains the metallic-type of conduction. Efros 15,21,43 
suggested that yes, all the c-liquid is metallic. Since the 
percolation through the c-liquid is now achieved in wide 
ranges of average v this suggestion lead him to the con- 
clusion that the peaks in <r xx must be wide, and cor- 
respondingly, the quantum Hall plateaux must be nar- 
row 43 (see also Refs. 42, 44). Indeed, this is what is 
observed at relatively high temperatures 45 . However, it 
is well-known that at lower temperatures peaks become 
narrow. One interpretation of these phenomena was sug- 
gested in Ref. 42 and is as follows. As the temperature 
goes down, some fraction of the c-liquid with filling fac- 
tors close to integer values becomes localized (or perhaps 
pinned); thus, the range of filling factors where c-liquid 
is metallic shrinks. If the want to preserve the property 
that only a single state per LLS is delocalized at zero 
temperature, then we have to assume that only very nar- 
row strips of the c-liquid where v(r) is half-integer remain 
metallic in the limit T —> 0. Therefore, we have to re- 
fine our criterion to have the conductivity peaks in the 
following way: the peaks corresponding to Landau level 
number N appear when there is a percolation through 
vn{v) = 1/2 (for spin up peak) and vn{v) = 3/2 (for 
spin down peak). This is a very strong assumption, but 
we do not know an alternative way to interpret the low- 
temperature magnetoresistance data. 

In our chess-board model the peaks in <r xx corresponds 
to the configurations of vn = 1/2 and vn = 3/2 fillings 
at the saddle points (Fig. 6). The difference in the av- 
erage filling factor for these configurations, which is bv, 
can be readily estimated as the ratio of the area occupied 
by the c-liquid near a single saddle point to the area of 
a chess-board cell, i.e., it is still given by Eq. (52) if by 
bx we now understand the size of the compressible area. 
Then as in the large density case, bv = 1/2 corresponds 
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to 6x ~ d £. Let us show that for such 6x the size of 
the incompressible region "1", which we denote by a (see 
Fig. 6a), is much smaller than 8x. 



a) 




b) 




FIG. 6. Fhe structure of a saddle point. Shown by grey 
is the area occupied by the c-liquid. Inside of this area 
two contours of constant filling factor, v(t) = 2N + 1/2 
(dashed line) and v(t) = 2N + 3/2 (dotted line), are im- 
portant for conductivity, (a) Average v is 2N + 1 and the 
same is v at the saddle point. Both contours are away from 
the saddle point. Fhus, they do not percolate through the 
sample and <r xx is zero, (b) Fhe average filling factor is 
decreased from the value of 27V + 1 by such an amount that 
v at the saddle point becomes 2N + 1/2. Fhis brings the 
former contour to the percolation condition, while the lat- 
ter contour is pushed further away from the saddle point. 
Fhis corresponds to the spin up peak in <r xx . 

We can estimate a in the way similar to the one used 
in Ref. 42 for the wire geometry. Without the exchange 
interaction the saddle point would be covered by a com- 
pressible region of size 6x and the external potential V(r) 
would be completely screened so that the electrostatic en- 
ergy be minimized. Due to the exchange interaction the 
region "1" appears. It can be considered as a quadrupole 
redistribution of the charge on top of the compressible 
region. This redistribution causes the deviation of the 
Hartree potential from its value for the complete screen- 
ing case. The magnitude of the deviation is of the order 
(e 2 riLa 3 /£ t )(i5j;)~ 2 I where the last factor is the typical 
value of the second derivative of the filling factor at the 
saddle point (for the complete screening case). Now a 
can be estimated from the condition that the deviation 



in the Hartree potential is of the same order of magni- 
tude as the deviation in the exchange potential, which is 
ahuj c if a £. Using Eq. (44) for the effective dielectric 
constant we find a ~ (i5a;) 2 / 3 (ai?) 1 / 3 . Replacing aR 
by £, we obtain 



2/3fl/3 



(54) 



Observe that for 6x £ we have £ <C a <C 6x, which 
makes our estimate consistent. 

For bv ~ 1/2 we have 6x ^ £ and, therefore, a <C 6x as 
we claimed above. In this case 6x as a function of N can 
be found assuming that the entire area of the saddle point 
is occupied by the c-liquid as it would be without the ex- 
change interaction. The evolution of the system at this 
stage is governed by the Hartree interaction (compare to 
Ref. 39), which makes the present case qualitatively dif- 
ferent from the large density case where the exchange 
part of the interaction played the major role. The reason 
for this difference is that in the present case the exchange 
interaction is no longer enhanced through the collective 
exchange field of neighboring saddle points — they are 
too far apart. 

Let us now calculate 6 function of N . Having 

done that, we will be able to find N c from the condition 
6x ~ d. Similar to the estimate of a earlier in this Sec- 
tion, that of 6x will be also based on the fact that the 
charge distribution in the compressible region can be de- 
scribed as a quadrupole. In this case it is the quadrupole 
formed on top of the "0"-"2" state. Indeed, in the present 
state there are more electrons in the two "0" regions, so 
there is a negative charge there, and there are fewer elec- 
trons in the two "2" regions, so they hold some additional 
positive charge. The size of the quadrupole is such that 
the sum of the Hartree and the external potential en- 
ergies has the lowest possible value. The energy of the 
quadrupole is, essentially, the product of its quadrupole 
moment and the second derivative of the external po- 
tential. Therefore, 6x can be estimated by equating the 
second derivatives of V and the Hartree potential: 



d 2 v 

dx 2 



e 2 n L 
s*6x ' 



which gives 



6x ~ jd, 



(55) 



(56) 



where j is the dimensionless parameter showing the rel- 
ative strengths of the interaction and the external poten- 
tial: 



1 



7 



e 2 n L 



d 2 V/dx 2 e*d 



(57) 



Eqs. (52,56) show that bv = 1/2 is achieved when j is 
of order unity. 

For the chess-board potential j is defined unambigu- 
ously because all the saddle points are identical. In the 
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random system this is not the case, but it is reasonable 
to assume that the point bv = 1/2 also corresponds to 
7 ~ 1 if in the definition of j we use a typical value of 
d 2 V/dx 2 , for which we take 



d 2 v 

dx 2 



2 1/2 

e a B n- 



d 3 



R 



1/2 



(58) 



[compare to Eq. (40)]. With the help of Eqs. (44,58) we 
can express j in terms of n, rii, and N: 



1/2 



3/2 



and hence, N c is given by 

N c ~ kydin/m) 1 ' 3 , 



(59) 



(60) 



which reproduces Eq. (4b) apart from the numerical co- 
efficient. The coefficient was derived in the following 
way. We performed a computer simulation for the chess- 
board potential given by Eqs. (50,51) with the purpose 
to find the ground state of charged liquid interacting via 
a Coulomb law with the dielectric constant . We im- 
posed an additional constraint that the density at every 
point is non-negative and does not exceed 2«l- The re- 
sults of the simulation depend only on two dimensionless 
parameters: j and v = n/riL- By varying the average 
density n in the system we could vary the density at the 
saddle points in the resulting ground state, in particular, 
to make it equal to and Sn^/2. We then calculated 

bv as the difference in v corresponding to these values 
of n. We found that bv = 1/2 is achieved at j = 0.1. 
Eventually, we substituted this number into Eq. (59) and 
obtained Eq. (4b). Note that the results for large and 
moderate densities [Eq. (4a) and Eq. (4b)] match at the 
crossover point defined by n/rii ~ a~ 3 , or more precisely, 
byEq. (5). 

Concluding this Section, we want to briefly discuss the 
case of small densities n/rii <C 1- This case is similar to 
the case of moderate densities with two exceptions. First, 
in this case R <C d and the averaged potential V(r) co- 
incides with the unaveraged one so that Eq. (40) has to 
be replaced by 



p(o)- 1 = £ 



2 1/2 

7r e z a B ?V 



and Eq. (58) by 



d 2 V e 2 a-Qn 



2„ I/ 2 



dx 2 d 3 
As a result, Eq. (56) gives 

fn\ 1/2 d 3 
bx ~ — — -. 
V rii J R 2 



(61) 



(62) 



(63) 



The other difference from the moderate density case is 
that the amplitude of V near the transition exceeds Tiu> c 



many times and v(r) varies typically not within 2N and 
2N + 2 but in a much wider interval. In this case the 
proper modification of Eq. (52) is 



bv 



(V 2 ) 1 ' 2 (b 



(64) 



Together with Eq. (63) this leads to Eq. (4c). The case 
of small densities is interesting in the sense that small 
regions (with size of order R) near the saddle points oc- 
cupied by the c-liquid determine, in fact, a large value 
of bv, which is a sample averaged quantity. So, bv can 
be of order unity while the average magnetization is very 
small. 

Our primary goal so far was to find N c where by defi- 
nition bv = 1/2. Another interesting question is the be- 
havior of bv for N > N c . For the moderate density case 
we can use Eqs. (52,56,59) to find that the transition has 
the tail in the form 



bv(N) ~ (N c /N) 3 , N >N C 



(65) 



Surprisingly, the same expression is valid for the small 
density case as well. In the derivation of Eq. (65) we 
assumed that bx > £. Now the question is whether or 
not it holds for smaller bx where the compressible region 
near the saddle point become more narrow than the re- 
gion "1" that it surrounds. To answer this question we 
have to include the exchange interaction in our consid- 
eration. However, Eq. (65) will be modified by at most 
a numerical factor (absent in Eq. (65) anyway) because 
at distances r <C £ the exchange interaction G ex (r), like 
the Hartree interaction, is equivalent to the Coulomb law 
with the dielectric constant [see Eqs. (15,44)]. 

Eventually, bx reaches /b- At this point our approx- 
imation of the system of discreet electrons by continu- 
ous liquid breaks down. For the chess-board bv should 
drop to zero. As for the random system system, there is 
a finite probability to find untypically soft saddle points 
(with small value of d 2 V / dx 2 ), which remain magnetized. 
Therefore, we can not exclude the possibility that bv is 
non-zero but is very small. 

Finally, we would like to remark that while for low- 
mobility samples or for high-mobility heterostructures 
with large electron densities it is possible to express N c 
in terms of the single-particle scattering time r, which 
is rather crude overall characteristic of the disorder, in 
the case studied in the present Section, N c is determined 
by the properties of some particular areas in the sample, 
namely, the saddle points of the disorder potential. 



V. CONCLUSION 



In this paper we demonstrated that the spin-splitting 
of the quantum Hall conductivity peaks is determined by 
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the competition between the disorder and the electron- 
electron interactions. We showed that in the mean-field 
approximation the distance bv between the spin-resolved 
<7 XX peaks measured in filling factor vanishes starting 
from a certain peak number N c if the bare Zeeman split- 
ting is zero. A non-zero Zeeman splitting smears the 
transition but in GaAs devices this splitting is very small 
and, in principle, can be totally eliminated by applying 
pressure. We calculated the peak number N c where the 
spin-splitting disappears as a function of the heterostruc- 
ture parameters and suggested a modified global phase 
diagram for the quantum Hall effect 25 , which now in- 
cludes spin (Fig. 2). Our calculations are in qualitative 
agreement with available experimental data 1,2,23 . How- 
ever, to verify our predictions in detail, experiments on 
high-mobility gated samples where the electron density 
can be varied are desirable 24 . Another way to verify these 
predictions would be the experiments where the effective 
density of randomly situated donors can be varied inde- 
pendently of the electron density (see Ref. 17). 

Our theory of the disorder-induced destruction of the 
many-body energy gap can be applied to other systems 
with additional degrees of freedom, examples of which 
can be multi-valley semiconductors or double quantum 
well structures. For instance, in Si MOSFET each branch 
of the phase diagram shown in Fig. 2 acquires an addi- 
tional "valley" fork (see also an experimental phase dia- 
gram by Kravchenko et a/ 46 ). 

The critical behavior, i.e., the question how bv goes 
to zero in the vicinity of the transition requires an ad- 
ditional investigation. Having in mind the experiments 
on gated samples where the electron density n can be 
varied while keeping the peak number constant, we can 
define the proximity to the transition by the dimension- 
less parameter (n — n c )/n c , n c being the density at the 
very transition. Our mean-field theory, which applies to 
the case of large electron densities and not too small bv, 
gives 

For moderate and small densities the dependence is more 
complicated (see the discussion at the end of Sec. IV). 
Note that the experimental study of the critical behavior 
may be difficult in view of the finite temperature effects. 
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